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ABSTRACT 

Aims. We aim at deriving the molecular abundances and temperatures of the hot molecular cores in the high-mass star-forming region 
NGC 63341 and consequently deriving their physical and astrochemical conditions. 

Methods. In the framework of the Herschel guaranteed time key program CHESS (Chemical Herschel Surveys of Star Forming 
Regions), NGC 63341 is investigated by using the Heterodyne Instrument for the Far-Infrared (HIFI) aboard the Herschel Space 
Observatory. A spectral line survey is carried out in the frequency range 480-1907 GHz, and further auxiliary interferometric data from 
the Submillimeter Array (SMA) in the 230 GHz band provide spatial information for disentangling the different physical components 
contributing to the HIFI spectrum. The spectral lines in the processed Herschel data are identified with the aid of former surveys and 
spectral line catalogs. The observed spectrum is then compared to a simulated synthetic spectrum, assuming local thermal equilibrium, 
and best fit parameters are derived using a model optimization package. 

Results. A total of 46 molecules are identified, with 31 isotopologues, resulting in about 4300 emission and absorption lines. High- 
energy levels (E u > 1000 K) of the dominant emitter methanol and vibrationally excited HCN (v 2 = 1) are detected. The number of 
unidentified lines remains low with 75, or < 2% of the lines detected. The modeling suggests that several spectral features need two 
or more components to be fitted properly. Other components could be assigned to cold foreground clouds or to outflows, most visible 
in the SiO and H 2 emission. A chemical variation between the two embedded hot cores is found, with more N-bearing molecules 
identified in SMA1 and O-bearing molecules in SMA2. 

Conclusions. Spectral line surveys give powerful insights into the study of the interstellar medium. Different molecules trace different 
physical conditions like the inner hot core, the envelope, the outflows or the cold foreground clouds. The derived molecular abundances 
provide further constraints for astrochemical models. 
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1. Introduction 

The astrochemical composition of the interstellar gas depends 
on the physical environment and the evolutionary state. In the 
cold dense regions molecules are unsaturated, while in con- 
trast saturated molecules are much more abundant in hot cores 
( |Herbst & van Dishoeck 2009). Unsaturated molecules are car- 
bon rich with long carbon chains, like radicals as C n H (n=2- 
8) or cyanopolyynes HC n N (n=3,5, 7,9,1 1), whereas saturated 
organic molecules are hydrogen rich with single bonds of car- 
bon like CH 3 OCH 3 . Molecular ions are predominantly unsat- 
urated in the cold phase, because most hydrogenation reac- 
tions are endothermic or hindered by potential barriers (Herbst 
1989 ). Ices on dust grains begin to form in the cold phase and 
are observed by broad absorption bands at mid-infrared wave- 
lengths. Sequential hydrogenitation of oxygen or nitrogen forms 
water or ammonia, respectively. Embedded in the water-ices 
are relevant amounts of CO, and through further hydrogena- 
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* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



dependent astrochemical models suggest that abundances of cer- 
tain molecules, such as H 2 S an d SO ? , might be used as an indi- 
cator of the evolutionary phase (Nomura & Millar 2004{ [Herpin| 
et al. 2009 ), and greater chemical complexity can be expected for 



longer evolution timescales (Garrod et al. 2008 1. Approximately 
165 molecules have been detected in the ISM or circumstellar 
shells. 

|Van Dishoeck ( 200 1 ) l emphasized the importance of spec- 
tral line surveys for the research of star-forming regions. Line 
surveys give the possibility to obtain a census of all atoms 
and molecules and give insights into their thermal excitation 
conditions and dynamics by studying line intensities and pro- 
files, which allows seprating different physical components. The 
discovery of new species or new excited lines from known 
molecules is possible. Furthermore, the wide range of energy 
levels and eventually opacities of one species allows a better 
determination of the temperature and abundance. The uniform 



1 A current list can be found at |http : / /www . astro . uni -koeln . | 
de/cdms /molecules 
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Table 1: Mass and temperature estimates for NGC 63341 



Source 


M(M Q ) 


Ti (K) p 


n (cm i ) 


Aral (Lq) 


Total" 


200 ± 100 


100 1.7 + 0.3 


1.2 x 10 7 


2.6 x 10 5 


SMA1* 


17 ± 50% 


100 






SMA2'' 


1 1 ± 50% 


100 






SMA3 6 


33 ± 50% 


60 






SMA4* 


36 ± 50% 


33 







Notes. Values taken from M |Sandell|j2000t and m |Hunter et aL]j2006f . 

M is the mass, Tj the average dust temperature, (3 the dust emissivity 
index, n the hydrogen density and L bo) the bolometric luminosity. 



recording with the same spectrometer in one observation reduces 
the calibration uncertainties. 

The giant molecular cloud NGC 6334, also known as the 
"Cat's Paw Nebula", lies in the Carina-Sagittarius Arm of the 
Milky Way at a distance of (1.7 + 0.3) kpc ( |Russeil et al.|2012| l. 
Studies at infrared wavelengths revealed sites of massive star 
formation (|Persi et al.|2005| ). NGC 63341 was first explored by 



Emerson et al. ( 1973[) asthe brightest infrared source in the 



northern part of the cloud, containing an OH maser and an ul- 
tracompact HII region. According to Emerson e t al.| ( fl973"] ), the 
pumping of the OH maser requires an energy source of a mas- 
sive collapsing protostar (M > 30 M G ). H II regions are associ- 
ated with high-mass stars (spectral classes O and B) where the 



UV radiation is strong enough to ionize hydrogen. Russeil et al. 
(2010) studied the densest cores in NGC 6334 in molecular line 
and dust continuum emission. For NGC 63341, they derive a total 
mass of M=1206 M G from the continuum flux density at 1.2 mm, 
an average density of n^-\2 x 10 6 crrT 3 and an upper limit for 
the virial mass of M v i r =525 M Q . The total mass estimate is rela- 
tively high, because a low value of the dust temperature of only 
20 K was assumed, as compared to the value of T=100 K derived 
by dSandell 2000), which leads to a factor of 6 difference of the 
mass estimate. While an average dust temperature of T=25 K is 
appropriate for most cores in the NGC 6334 complex (Matthews 



et al.|2008] l, the dust temperature is known to be much higher in 
NGC 63341. 

Interferometric studies revealed that NGC 63341 consists of 
several cores, labeled SMA1 to SMA4 ^Hunter et al.|2006| >. The 
separation distance between SMA1 and SMA2 is about 3.5" or 
6000 AU. SMA1 and SMA2 show rich spectra of molecular tran- 
sitions and SMA3 coincides with the H II region which exhibits 
free-free and dust emission and is excited by the near-infrared 
source IRS IE ( |De Buizer et al.|2002[ ). SMA4 only shows dust 
emission but no line emission, which is still unexplained. Table[T] 
summarizes characteristic values for NGC 63341 derived from 
the dust continuum emission at different wavelengths and, which 
classifiy the source as a low-luminosity high-mass hot core. The 
overall size of NGC 63341 is 10" x 8", where 1" equals to 
1700 AU. SMA1 and SMA2 are so-called hot molecular cor es, 
where certain attributes apply according to Cesaroni (2003}: a 
high temperature (T >100 K), sizes smaller than 0.1 pc, large 
masses (10 - 1000 M Q ) and luminosities exceeding 10 4 L Q . 
Furthermore they are often associated with water masers and ul- 
tracompact H II regions. Hot cores are heated by high-mass pro- 
tostars which are embedded in a dusty gas envelope. Their rich 
molecular spectra originate from the evaporation of ice mantles 
covering the dust grains, which release at T ~ 100 K molecules to 
the gas phase. Molecules released from grain mantles then begin 
to drive a fast, high-temperature gas-phase chemistry, forming 
complex organic species like acetone, beginning from precursors 
like methanol. Their rotational, vibrational and radiative excita- 




Right Ascension (J2000) 

Fig. 1: Integrated SMA line intensity map over all species. 
Overlaid in black contours is the continuum emission. The con- 
tour levels are from 90% to 10% of the maximum in steps of 
10%. The large and small white circles represent the FWHM 
Herschel beam size for band 1 and band 7, respectively. 



tion leads to a characteristic emission spectrum. This explains 
why the spectrum is poorer in earlier stages like the colder pro- 
tostel lar core NGC 6334I(N) ( |Brogan et aLl|2"0"0"9{ |Walsh et al. 
2010]), 2' or 1 pc north of NGC 63341. 



Multiple spectral line surveys have been carried out for 
NGC 63341 : |McCutcheon"eT al. (2000) 334-3 48 GHz, |Schilke 
|eTaLl(|2006l ) 459-461 GHz and 817-81 9 GHz, |Thorwirfh et al. 
( 120031 ) 88-115 GHz and 218 -267 GHz,|Walsh et al.|d20101 i 84- 
116 GHz and the latest from Kalinina et al. (2010) in the range 
81-242 GHz. Besides, several maser lines have been observed 



from NH 3 (Beuther^TaT1[2007]l, CH 3 OH (Walsh et al.||1998i 



and H 2 (Migenes et al. 1999) as well as bipolar outflows in 
CO ([Leurini et al.||2006[ |Qiu et al.||20T"fl ) and SiO ( [Gibb et al.| 
2007). This source is part of the CHESS program (Chemical 
Herschel Surveys of Star Forming Regions), a key program of 
the Herschel Space Observatory ( |Ceccar elli et al. 2010 ). The aim 
is to study eight different sites and phases of star formation vary- 
ing in mass, luminosity, evolutionary state, astrochemical com- 
position etc. by conducting and comparing their spectra. In this 
article, the entire spectrum (500-1900 GHz) of NGC 63341 is 
analyzed which was observed by Herschel in 2010. Previously, 
some small sections of the spectrum were studied to find tran- 
sitions of new molecules not observed before in the interstellar 



medium (ISM): Oxidanium yl (H,Q + ) by Ossenkopf et al. (2010) 
and chloronium (H 2 C1 + ) by| Lis et aL|(|2010|. A de tailed analysis 
of water has been made by Emprechti nger et ak] ( |2010| ) and of 



methylidyne (CH) by | van der Wiel et al.| ( |2010| ). 



One difficulty is that if a large beam encompasses a complex, 
unresolved source structure, the spectrum can become complex 
due to different overlapping components. The spatial informa- 
tion from the SMA is complementary. The Herschel beam is 
larger than NGC 63341 so that the source is spatially unresolved, 
while the SMA beam can resolve the individual cores. The SMA 
maps thus reveal the source morphology and give insights into 
the distribution of various molecules. This helps to decompose 
the contributions from different components present in the HIFI 
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Fig. 2: Complete HIFI spectrum of NGC 63341. B-bands are coloured in blue and bands 5-7 are smoothed by averaging over 5-10 
channels. 



spectrum and to differentiate between the two hot cores SMA1 
and SMA2, the extended emission and outflows, see Fig. [T] 



2. Observations 

2.1. Observations with Herschel 

Observations were conducted with the high spectral resolu- 
tion receiver system HIFI (Heterodyne Instrument for the Far- 
Infrared; Ide Graauw et al. 



Table 2: Overview over all HIFI bands for the observation of 
NGC 63341 



Observatory (|Pilbratt et alT 



2010) aboard the Herschel Space 



2010 1. HIFI is divided into seven 



bands with different receivers which overlap at the band edges. 
Each band is separated again into two parts a and b us- 
ing different polarizations and local oscillator chains. The 
full frequency range is from 480-1907 GHz, with a gap be- 
tween 1280 and 1410 GHz. The first five bands employ SIS 
(superconductor-insulator-superconductor) mixers with an inter- 
mediate frequency (IF) bandwidth of 4 GHz. Bands 6 and 7 are 
based on Hot-Electron Bolometer (HEB) mixers and have an 
IF bandwidth of 2.4 GHz. The HPBW of the 3.5 m Herschel 
telescope varies with wavelength from 41" to 12". Table 
gives a log of the HIFI observations of NGC 63341 and Fig. 
shows the whole calibrated spectrum. The continuum level in 
creases continuously with frequency, except for band 5b, where 
a small offset is seen. The line emission in band la contributes 
only 3% to the total flux because the HPBW at 520 GHz is 
four times bigger than the source size so that in addition the 
continuum from the surrounding area is observed. The emis- 
sion lines with the highest peak intensity, T >15 K, are those 
of CO and 13 CO. The telescope was pointed toward the coor- 
dinates o-(J2000): 17 h 20 m 53.32 s , c5(J2000): -35°46'58.5", be- 
tween SMA1 and SMA2. 

The observations were carried out mainly in two sessions 
in March and October 2010, in dual beam switch mode with 
an internal chopper mirror switching between the source and 



Band 


Freq. range 


Beam size 


Date 


Obs. time 


rms 




(GHz) 


(") 




(h) 


(K) 


la 


480-560 


41 


28.02.2010 


3.6 


0.015 


lb 


554-636 


02.03.2010 


2.2 


0.03 


2a 


626-726 


29 


09.10.2010 


3.7 


0.04 


2b 


714-801 


19.03.2010 


3.3 


0.05 


3 a 


799-860 


25 


09.10.2010 


1.8 


0.05 


3b 


858-961 


19.03.2010 


3.9 


0.07 


4a 


949-1061 


21 


09.10.2010 


4.8 


0.09 


4b 


1052-1122 


03.03.2010 


2.9 


0.10 


5 a 


1108-1244 


18.5 


13.10.2010 


4.2 


0.37 


5b 


1227-1280 


12.03.2011 


1.7 


0.37 


6a 


1410-1575 


14.5 


14.10.2010 


0.9 




6b 


1575-1700 


14.10.2010 


3.6 


1.20 


7a 


1703-1799 


12.5 


14.10.2010 


2.6 


1.27 


7b 


1788-1907 


13.10.2010 


3.2 


1.27 



two reference positions 3' away from the target to provide flat 
baselines and calibrate the continuum. The total observing time 
was ~42 h. The wide band spectrometer provided a spectral res- 
olution of 1.1 MHz, resulting in a velocity resolution of 0.7- 
0.2 kms . The double sideband spectra were processed with 
HIPE (Herschel Interactive Processing Environment; Ott 2010) 
version 5, deconvolved to single-sideband spectra following the 
algorithms described in ( {Comito & Schilke|2002 1, and exported 



to CLASS (see Sect. 3. 1 1 for subsequent analysis. All spectra are 
corrected for the beam efficiency and the intensity scale is the 
main-beam brightness temperature in K. The spectra presented 
here are equally weighted averages of the horizontal and verti- 
cal polarizations, in order to reduce the noise. Band 4a has some 
sinusoidal baseline features from standing waves. For band 6a 
only a small portion around a HC1 + transition was observed, but 
the data are corrupted and not discussed here. The band cover- 
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age is unbiased, meaning that all frequency ranges are included 
with the band edges where the noise increases. Different side- 
band gains can lead to discrepancies between line intensities in 
the two bands in the overlap region, see Fig. |B.lg1 for an exam- 
ple in the J = 6 - 5 transition of N 2 H + at 558.966 GHz. The 
discrepancies amount to 10-20% in the bands 1-3 and become 
negligible in the higher bands. In the further analysis, this un- 
certainty can be avoided and has marginal effects when several 
lines of a given molecule are present. 

2.2. Auxiliary observations from the Submillimeter Array 

Supporting SMA data at 230 GHz are kindly provided by 
C. Brogan and T. Hunter. The continuum map is identical to 
that published in |Hunter et al.| ( [2006[ ), whereas the spectral im- 
age cubes are unpublished. Detailed information about the ob- 
servation and calibration can be found therein. The data were 



calibrated again in 2007 in Miriad (described in Brogan et al. 
2009). The separation of the line and continuum emission was 
performed in the uv plane. For the purpose of modelling, the 
interferometric line and continuum emission were recombined 
in CLASS for the portion of the spectrum containing the fore- 
ground CN absorption feature. The observed frequency range 
covers two 2 GHz wide bands, 216.6-218.6 GHz in the lower 
sideband and 226.6-228.6 GHz in the upper sideband, with a 
frequency resolution of ~0.8 MHz. The map covers a 64" x 64" 
field and the synthesized beam is 2.2" x 1.3" in the USB. The 
line-to-continuum ratio is 12% for SMAl and 21% for SMA2 in 
the given bands. 



3. Analysis: line identification and modeling 

3.1. Simulating spectra with XCLASS 

A spectral line is identified by comparing the center frequency 
with a measured frequency of a known chemical species (atom, 
molecule, ion) in the laboratory. Two large databases relevant for 
molecular astrophysics, from the radio to the far-infrared region, 
have been built up in the last decades: The Cologne Database for 
Molecular Spectroscopy (CDMS) ( jMiiller et al]|2005[) and th e 
Jet Propulsion Laboratory (JPL) catalogue ( |Pickett et al.|1 998). 
If the exact source velocity is unknown or several candidates are 
present for the identification of a spectral line, secondary criteria 
such as the upper energy level or the Einstein A coefficient must 
be taken into account. Although many lines in the spectra of hot 
cores are blended, the identification in spectral surveys is more 
confident because more lines of a species are present, which can 
then be checked for consistency, i.e. all lines of a molecule which 
are expected to be there should be found with the expected line 
strength, line shape and velocity. 

The main program used in this work is XCLASS, which was 
described in |Comito et al. ( 2005 \ and applied to a ground-based 
Orion KL line survey. It is an extension of CLASS (Continuum 
and Line Analysis Single-dish Software), which is part of the 
GILDAS0 package. XCLASS includes entries from the CDMS, 
JPL and a private catalogue, which contain all necessary infor- 
mation: rest frequency v/, integrated intensity, lower state energy 
Ei, upper state degeneracy g up , quantum numbers and the parti- 
tion function Q. With this database, it is possible to simulate a 
spectrum of a molecule, by providing parameters such as the ex- 
citation temperature and abundance, and compare the simulated 
spectrum with the observed one. Similar software tools devel- 



oped in the last years are Weeds ( Maret et al.|2011| l, a package 
for CLASS, and CASSIS ( |Caux et al.|201 1) . The hne identifica- 
tion for NGC 63341, which transfers to a molecule identification 
in XCLASS, is made easier because 20 detected molecules are 



reported in the line survey from |Kalinina et aL ( 2010| l at lower 
frequencies and in other surveys mentioned in Sect.[l] 

The function which models a spectrum in XCLASS is a solu- 
tion of the radiative transfer equation for an isothermal object in 
one dimension, called detection equation ( Stahler & Palla 2005} . 
Furthermore, the finite source size and dust attenuation are con- 
sidered. The model function is: 



T(v) 



XiZ ??( ^ c) f- /(r " 



c )-J(T bg )](l-e- T ^""')-e- 



(1) 



with the optical depth 



T(vr c = ^ T ( v > 



I 



f \tn,c 

t(v) ; ' = 



8ttv 3 



A U ,N" 



Q(m, T ex ' ) 



(1-- 



(2) 

-Av/*T« ) ( y ) m )C , (3) 



the dust optical depth 



Td= ^'^-(2^Hzf- 33Xir26 ' 



the Gaussian line profile 



and the beam filling factor 

rj(fim,c) 



exp 



( (v-(v/ + v LSR )r x 

2<T 2 



(4) 



(5) 



(6) 



http : //www . ir am . f r/IRAMFR/GILDAS 



The radiation temperature J(T) is defined as J(T) = j-j^/wzf- 
The sum goes over the indices m for molecule, c for compo- 
nent and / for spectral line transition. 7j„ is the background 
continuum temperature and 9 t the telescope beam size. The line 

("CO 

profile function is normalized so that J dvtp(v) = 1. The fac- 
tor e~ Td takes the dust attenuation into account where the dust 
mass opacity k is taken from Ossenkopf & Henning (1994). 
For NGC 63341, we adopt a value of /? = 2 and Ki.3 mm = 
0.42 cm 2 g _1 , which is between the values for non-coagulated 
dust grains with thin ice mantles (/ci.3 mm = 0.51 cm 2 g _1 ) and 
without ones (k\ 3 mm = 0.31 cm 2 g _1 ). The total hydrogen col- 
umn density N H = A^(H)+2A^(H 2 ) is 3 x 10 24 cm -2 , which is 
derived from the average hydrogen density given in Table [T] and 
a source size of 10" and is in agreement with values from previ- 
ous articles ( |van der Wiel et al.|2010|[Lis et al. 2010 1 considering 
the uncertainties. 

Overall, there are five input parameters for one component 
of a molecule: 

m , c source size in arcsec 

T ex excitation temperature in K 

iV to t total column density in crrr 2 

Av velocity width (FWHM) in km s -1 

v L sr source velocity in km s _1 . 

Components can be identified as spatially distinct sources within 
the beam, such as clumps, hot dense cores, colder envelopes or 
outflows, and can usually be distinguished by different radial ve- 
locities. They do not interact with each other radiatively, but are 
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simply superposed in the model. However, for absorption lines 
the emission from other components is considered first by cal- 
culating the background emission spectrum, which is then used 
as a new "continuum" for absorption lines. By fitting all species 
and their components at once, line blending and optical depth ef- 
fects are taken into account. Because of the high densities in the 
hot cores (n^ > 10 7 cirT 3 ), it is assumed that most molecules 
are thermalized at the gas temperature, so that the assumption 
of LTE is reasonable. With the LTE condition, every popula- 
tion of a level i is known and given by the Boltzmann distri- 



bution «,■ = %gie 



-E,/(kT) 



The modeling of a spectrum speeds up 



line identification enormously, since several hundred lines of a 
molecule like methanol are predicted at once in a HIFI band. 

The modeling can be done together with isotopologues 
(and higher vibrational states) of a molecule assuming a 
conversion factor stored in a separate database. All param- 
eters are expected to be the same except the abundance. 



For carbon, [Milam et al. (2005 ) derive an average gradient 
of 12 C/ 13 C=6.21(1.00) • D cc + 18.71(7.37), where D cc is 
the distance to the Galactic center in kpc. The carbon iso- 
topic ratio for NGC 63341 would be 57(10). Eight isotopic 
species are found in this survey, besides the main iso- 
topes. The ratios used for NGC 63341 are: [ 35 C1]/[ 37 C1]=3, 
[ 32 S]/[ 34 S]=23, [ 32 S]/[ 33 S] = 127, [ 12 C]/[ 13 C]=60, 

[ 14 N]/[ 15 N]=300, [ I6 O]/[ I8 O]=500 and [ 16 O]/[ 17 O] = 1500. 
Deuterium is found also, but here no common ratio is taken 
because the chemical fractionation is very high. The ortho to 
para ratio for some molecules is assumed to be the statistical 
equilibrium value (high temperature limit). 

3.2. Model fitting with MAGIX 

Once an initial guess is set up, the parameters can be fitted using 
the program MAGIX (Modelling and Analysis Generic Interface 
for eXternal numerical codes; |Bernst et al.||201 lj |Moller et al.| 
2012| l. It delivers a package of different optimization algorithms 
and an interface to XCLASS. For every of the five fit parame- 
ters a range can be set over which to search for the best fit, or 
the parameter can be kept constant. Used here is the Levenberg- 
Marquardt (LM) algorithm with the method of least squares, 
which tries to minimize the squares of the difference between 
the fit function y(Xj\ao . . Mm-i) with the M fit parameters a and 
the data. The sum over all n data points (x^yi) is defined as: 



y 0; - y(Xj\a . . . a M -l)) 2 



err 



(7) 



The uncertainty for y, is cr,-. The iteration stops when a min- 
imum is found and the algorithm has converged, or when the 
change in x 1 from one iteration step to another is very small, 
e.g. Ay 2 < 10~ 6 . Since this function is a priori unknown, it be- 
comes clear that the algorithm can converge into a local and not 
a global minimum. A bad fit for all parameter values can be a 
hint for non-LTE populations, like HCN lines with high critical 
densities, or maser activity. The automated fitting has to be cor- 
rected manually sometimes, since different problems can occur, 
which are described in the following. 

Starting values: The choice of the initial values for the pa- 
rameters, is fundamentally important in the LM algorithm. Bad 
or random starting values lead to an unnecessary amount of com- 
putation time. The starting values are set manually by the user 
and should be physically reasonable. Because of their fit be- 
havior, the five fit parameters in XCLASS are divided into two 
groups: 6,T ex ,N and Av, vlsr- The last two are generally well 



fixed and do not change anymore, so they are held constant. For 
high opacities, the line width is broadened which is known as op- 
tical depth broadening. To check for a local or global minimum, 
MAGIX was iterated again with another set of starting values 
and then checked if it still converged to the same point. The dis- 
advantage is that this takes computation time with several repe- 
titions and the variation possibilities increase dramatically with 
more velocity components for a given molecule. This check has 
been done regularly with molecules with one component, but in- 
frequently for more complicated cases. The best fit parameters 
for the species with multiple overlapping components may thus 
not be uniquely determined. 

Calculation time: The time increases with the number of 
parameters and data points. To reduce the calculation time, es- 
pecially for "weeds", not all transitions from every band are 
included but a variety of optical thin and thick lines, lowest 
and highest frequencies and accordingly a big range of differ- 
ent energy states. Suitable for this are Q-branches of methanol 
or dimethyl ether. Weeds is the colloquial term for interstellar 
molecules, which show a plethora of spectral lines. They are all 
organic asymmetric rotors with a methyl group (-CH 3 ), except of 
S0 2 . The classification is not strict, but these five are considered 
to be the most important ones: methanol, dimethyl ether, ethanol, 
methyl formate and sulfur dioxide. In methanol, CH 3 OH, the 
methyl group can tunnel the potential barrier at every 120° rel- 
atively to the hydroxyl group (-OH), leading to an internal rota- 
tion or torsion about the C-0 bond. The torsional levels are split 
into sub-levels (A and E symmetry), and the interaction of tor- 
sional and rotational motions yields multiple energy states and 
a complex spectrum. Furthermore, heavier and larger molecules 
have a larger moment of inertia, which is inverse to the rota- 
tional constants, so that they produce denser spectra. The deter- 
mination of the excitation temperature with XCLASS is much 
faster than the so called Boltzmann plot or rotational diagram 
method, where the upper state abundance N u of a line transition 
can be derived from the integrated area in the optical thin limit. 
This method has the disadvantages that it is unpractical for many 
lines and misleading for line blends and high optical depths. 

Uniqueness: The model function Eq. [T]has two limits where 
degeneracy occurs: in the optically thick limit t s> 1, the last 
term vanishes and J{T ex ) and rj(6) enter linearly into the equa- 
tion; in the optically thin limit r <sc 1, (1 - e~ T ) ~ t and rj(9) 
and N enter linearly. This results in the fact that these two pa- 
rameter pairs cannot be determined independently. In the fitting 
procedure this can be noticed in the behavior that x 2 changes 
only slightly while one of the parameter increases and the other 
decreases. It is not always obvious if a spectral line is optically 
thick, especially when several components add together. A solu- 
tion to the problem of degeneracy are isotopologues, if present. 
With the ratios mentioned above, a consistent model can be fitted 
with constraints on N and also for Av. In some cases the rare iso- 
topologues are fitted well, but the main isotopologue is optically 
thick and too weak to explain the high observed intensities; this 
is the case for CH 3 OH, SO, H 2 CO and HCO + . This can be a hint 
that a second component is present when no signs for that could 
be deduced from the line shape. However, every additional emis- 
sion component should be justified, since this makes the model 
more complex and increases the possibility of redundancy. In the 
case of NGC 63341, two near overlapping components are justi- 
fied by the SMA observations which show that both cores SMA1 
and SMA2 show line emission at different velocities. 

Weighting: x 2 is the sum over a U modeled spectral lines. 
One weighting could be done by the uncertainty of y,, <r,, which 
is currently not implemented in MAGIX. But this is not a severe 
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Fig. 3: (a) Fit example around 497 GHz showing the data (black) and the fit (red). The decomposition into the contributions of 
individual species is shown above in (b). 



limitation, because the noise or rms is nearly constant across a 
given HIFI band. It is relevant when the line intensity of the main 
isotopologue is much higher than for the less abundant isotopo- 
logue. Stronger lines are then more weighted compared to weak 
lines because the change in^ 2 is larger. This can lead to a model 
that is biased toward fitting the stronger lines. Another problem 
is the fitting in higher-frequency bands with larger noise (band 
5-7). Here the algorithm tends to increase the line width because 
it smooths the line wings better. Line widths from transitions in 
lower bands reveal then that this is an overestimation. To fix this, 
the line width must be held constant and the best fit value is taken 
from transitions in lower bands. 

Error estimate and confidence limits: It is assumed that 
in general the standard deviation of 6, T and N is at least 10%, 
because a small higher or lower value does not change the fit sig- 
nificantly. Additionally, there are other errors not accounted for 
in the fit like calibration uncertainties, uncertain isotopic ratios, 
N(H2), etc., so that in the worst case only the magnitude of N is 
roughly known. 



Due to the fact that lines of one species are blended with lines 
of other molecules, the fitting procedure is done in three major 
steps: 

1. Fitting intense lines and weeds. By modeling weeds one 
filters them out and takes into account the line blendings. For 
intense lines of diatomic molecules like SO, the effect of line 
blending is small. Blended lines are first excluded from the fit- 
ting range. 

2. Fitting weaker and absorption lines. For weak lines, the 
noise becomes important so that the fitting procedure is less re- 
liable. Too weak lines, below 2cr, of a molecule are therefore 
excluded from the fitting range, except for cases where they can 
serve to give an upper limit for T or N. Absorption lines are not 
expected to be associated with the hot cores in NGC 63341 be- 



cause they come mainly from cold foreground gas. Apart from 
that are ground state transitions with low energy levels from H 2 
or NH 3 . 

3. Repeating all in a second iteration. The difference now is 
that the simulated spectrum of all molecules are used as a new 
background. The parameters are saved globally in a file and read 
in by MAGIX, but only the molecule of interest is fitted again. 
This time, blended lines of a species are taken into account and 
more components added, if necessary. A simultaneous fit of all 
parameters of all components of all molecules is in principle pos- 
sible, but far too time consuming and inefficient. In the end, the 
result is an improved and consistent fit, see the extract in Fig. [3 
This could be done iteratively several times, but it was noticec 
that the changes are marginal and larger only for molecules with 
very few transitions present. 

3.3. Analysis of interferometric data 

The presence of several components with small radial velocity 
differences which overlap makes the assignment to one certain 
physical source difficult. The advantage of interferometer data 
is having the spectral and spatial information together. On the 
other hand, the SMA frequency range of 4 GHz is small com- 
pared to HIFI and extended emission can be missed, because it 
is filtered out by the interferometer. Two spectra were extracted, 
one from the peak emission toward SMA1 and the other toward 
the peak emission of SMA2. It must be noticed that for some 
molecules like methanol, the emission overlaps from both cores. 
A direct comparison between them reveals that firstly SMA1 
is constantly red-shifted compared to SMA2, secondly the line 
width in SMA1 is broader and thirdly some lines are only present 
in SMA1 or SMA2. The intensities in Jy/beam are converted to 
Kelvin by using the formula 

1.22 x 10 6 -F(Jy) 



7mb(K) = 



[v(GHz) • 0(")] 2 



(8) 
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where 6 is the beam size. The conversion factor is then about 
12.5 K/Jy for the LSB and 8.5 K/Jy for the USB. The rms is 
about 2.5 K in the LSB and 2 K in the USB. Maps of the inte- 
grated intensity over one spectral line of all detected molecules 
were made with the toolkit Karma ( |Gooch|1 996) and the source 
sizes were derived from Gaussian fits. The whole procedure ex- 
plained in Sect. |3.2| was performed again with the SMA spectra. 
The model derived from the HIFI spectrum was extrapolated to 
the lower frequencies of the SMA bands and compared with the 
observed one. Deviations can be explained by the fact that the 
HIFI model includes the outer parts like the envelope and is an 
average over a range of temperatures, whereas the SMA spectra 
correspond to the inner hot core with higher temperatures. 

The new derived source sizes 6 were adapted for the HIFI 
spectrum. For some molecules, the source sizes derived from 
the HIFI spectrum were already in good agreement with the one 
from the SMA data. Small deviations (< 2") are permitted when 
they lead to a lower x 1 value. The source size is likely to change 
when excitation conditions differ, so that the high-/ transitions 
in the HIFI spectra may come from a more compact region than 
that observed in the SMA maps. For the cases where the lines 
are optically thin or thick, 9 is held fixed in order to cancel out 
the ambiguity. 



4. Results 

4. 1 . General results 

A total of 41 molecules with their isotopologues are identified 
in the HIFI spectrum. The best fit results are given in Table [3] 
The whole survey would need several hundred pages of fig- 
ures and tables with line lists. Therefore, to give an impression, 
only selected fits for each molecule are attached in Appendix B. 
The modeled spectrum for each molecule and the full observed 
and modeled spectrum for NGC 63341 will be made available 
online, after reprocessing the HIFI spectra with HIPE 8. See 
Appendix B for further details. The average source velocity 
(considering all components in the range -5 to -10 kms 1 as- 
sociated with NGC 63341) is -7.3 km s~' and average line width 
4.2 kms -1 with a standard deviation of about + 1 kms -1 for 
both spectral features. The velocity width is dominated by tur- 
bulence, since the thermal Doppler width contributes less than 
1 kms -1 for typical values of 7=100 K and m=20 amu. Some 
molecules show strong broad wings which are produced by out- 
flows: CO, CS, SiO, H 2 and HCN. H 2 CS has an asymmetric 
line profile with a stronger redshifted wing. Especially the com- 
ponent at -17 to -20 kms -1 , observed in C, NH 3 and H 2 S with 
T - 10 - 30 K is presumably the blue outflow lobe. Two cold 
(~5 K), extended foreground clouds are detected consistently 
with components around -2 km s~' and +7 km s" 1 in the follow- 
ing species: C and C + , CH and CH + , H 2 and H 2 + , HF, H 2 C1 + 
and OH + . This is in agreement with an observed OH absorp- 
tion component at +6 km s^ 1 towards several NGC 6334 cores 
reported by |Brooks & Whiteoak1 ( |2001| l. 

The best fit results for the SMA spectra are presented in 
Table A.l the list of detected lines in Table A. 2 and the spec- 
tra for the LSB and USB in Fig. |A.l| in Appendix A. Overall, 20 
molecules are identified, of which 15 are also present in the HIFI 
observations. The additional five molecules are all organic: ace- 
tone (CH3CH3CO), ethylene glycol ((CH 2 OH) 2 )), formic acid 
(HCOOH), acrylonitrile (C 2 H 3 CN) and propionitrile (C 2 H 5 CN). 
The latter two only have low intensity features (< 3<x) in the 
HIFI band la spectrum and are therefore considered to be not de- 
tected. Also, 3 isotopologues of OCS are found. In general, the 



fits give higher temperatures than those derived from the HIFI 
data, up to 200-300 K, because the emission originates from the 
center of the hot cores, observed in a smaller beam. Although 
not true for all molecules, SMA2 has on average higher ex- 
citation temperatures, lower column densities, abundances and 
smaller line widths than SMA1. On average, the source velocity 
forSMAl is -5.6i0.lkms- 1 andforSMA2 -7.6 + 0.1 kms" 1 , 
so that the relative radial velocity difference is 2 kms -1 . When 
applied e.g. to ethanol (C 2 H 5 OH), which is only seen toward 
SMA2, the source velocity derived from the SMA spectra is con- 
sistent with that from HIFI. In the case of C 2 H 4 or NH 2 CHO, a 
small discrepancy of 0.8 kms 1 is seen, which can be explained 
by the fact that in the SMA bands only a few transitions are de- 
tected, some of which are blended, so that the precise velocity 
determination is hindered. 

The chemical differences between both cores are deter- 
mined from the SMA maps and the velocity components in the 
HIFI spectra. Those molecules that show dominant emission or 
only emission in SMA2 are HCOOH, (CH 2 OH) 2 ), NH 2 CHO, 
C 2 H 4 and C 2 H 5 OH. For SMA1, the dominant species are 
HC 3 N, C 2 H 3 CN, C 2 H 5 CN, H 2 S, CN, HCN and HNC. It is 
noticeable that more oxygen-bearing species, including water 
(Emprechtinge r~et al.|2012| l, arise in SMA2 and more nitrogen- 
bearing species in SMA1. This variation in the peak abundance 
is also observed in other star-forming regions, such as Orion KL 
between the hot core and the compact ridge (Comito et al. 2005 ). 



4.2. Statistics 

A spectral line is considered to be detected when the peak in- 
tensity exceeds the 3<x level of the random noise. Furthermore, 
the line width should be comparable to that of all other detected 
lines. In this way one excludes single channel spikes and sets a 
lower limit for the integrated line intensity. A lower limit of 2cr 
would be acceptable when this limit is not taken for a single line 
of a molecule but rather for many lines since the probability that 
this is a coincidence decreases when having many lines observed 
and a good model fit. Nevertheless, some weak features are in- 
cluded in certain cases in the fitting range because they can give 
useful upper limits for the model. 

Overall, 49 chemical species have been found in NGC 63341 
so far, listed in Table|4] ranging from fine structure transitions of 
C, through over mostly di- and triatomic molecules, to organic 
molecules with 10 atoms. About half of them were found in pre- 
vious line surveys. Methyl acetylene (CH 3 CCH) has been iden- 
tified by (Kali nina et al.|2010[ ), but is too weak in HIFI (< lcr). 
We find 3 1 isotopologues in our data. The number of lines in all 
HIFI bands derived from the XCLASS fits is 4217, see Table|5] 
Hyperfine structure splitting is not considered and counted as 
one line. For the number of unidentified lines see next section. 
An automated line counting on the observational spectra gives 
the same number in orders of magnitude, but is too inaccurate at 
the band edges and for blended lines. The integrated line inten- 
sities are derived from the simulated spectrum, and for bad fits 
like CO, CS, HCN etc. from Gaussian fits. Because the fits are 
not perfect, the presented values should be viewed with caution. 

The line density is 5.5/GHz in the first 5 HIFI bands, ranging 
from 17/GHz in la to <l/GHz in 5a. Bands 6 and 7 have a high 
rms and less than 10 lines each. It is evident that the five weeds 
mentioned in Sect. |3.2| account for over 75 % of all lines, and that 
methanol is the dominant emitter besides CO. Methyl formate 
has a very dense spectrum, but is less intense at HIFI frequencies 
than in both SMA bands where it has the highest number of lines. 
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Table 3: Fit results for the HIFI spectrum of NGC 63341 
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Notes. Several rows for one species denote multiple components. Extended source sizes are abbreviated with ext. (1) including 13 CH + . (2) including 
15 NH 3 . (3) only isotopologues H 2 s O, H 2 7 0. (4) only isotopologues 13 CO, C 18 0, 13 C ls O, C 17 0. (5) only isotopologues H 13 CN, HC 15 N, H 13 C 15 N. 
(6) including H 37 C1. (7) including H 3 , 7 Cl + . (8) including H? 4 S. (9) only isotopologues C 34 S, C 33 S, 13 CS, 13 C 34 S. (10) including 34 SO, 33 SO. (11) 
including H 2 C 34 S. (12) including l3 CH 3 OH. (13) including 13 CH 3 OCH 3 . (14) including H 2 3 CO. (15) including H 13 CO + . (a) t small, (b) outflow 
component. 



Table 4: List of all detected chemical species in NGC 63341 
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4.3. Unidentified features 

Table|6]lists the frequencies and strengths of the 75 unidentified 
lines in the HIFI spectra and the 27 unidentified lines in the SMA 
spectra. The strength of the lines are derived from Gaussian fits. 



211 unidentified features were detected initially, going up from 
band la to 4b, and none in the HEB Bands. All of them are emis- 
sion lines except of one weak absorption line at 507 686 MHz 
below 5cr. Most lines are relatively weak. The highest peak-to- 
noise ratio of 14 has the line at 492 325 MHz. Erroneous features 
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Table 6: List of unidentified lines in the SMA and HIFI spectra for NGC 63341 



Frequency 


Area 


I 


I/rms 


Width 


Frequency 


Area 


I 


I/rms 


Width 


(MHz) 


(Kkms- 1 ) 


(K) 




(kms-») 


(MHz) 


(Kkms" 1 ) 


(K) 




(km s" 1 ) 


SMA 




216671 


41.2 


8.1 


3.2 


4.8 


227128 


68.0 


11.3 


5.7 


5.7 


216715 


71.7* 


13.5 


5.4 


5.0 


227239 


26.0* 


8.1 


4.1 


3.0 


216913 


48.0 


12.4 


5.0 


3.6 


227244 


32.3* 


7.0 


3.5 


4.3 


216919 


82.8 


24.1 


9.6 


3.2 


227534 


61.0 


14.2 


7.1 


4.0 


217000 


54.9 


11.4 


4.6 


4.5 


227728 


43.7 


12.8 


6.4 


3.2 


217056 


90.0 


20.8 


8.3 


4.1 


227803 


41.6 


11.4 


5.7 


3.4 


217262 


91.5* 


22.3 


8.9 


3.9 


227823 


47.8* 


7.4 


3.7 


6.1 


217427 


38.2 


11.6 


4.6 


3.1 


227841 


99.0* 


21.1 


10.6 


4.4 


217496 


70.9 


18.8 


7.5 


3.5 


228001 


53.0 


7.6 


3.8 


6.6 


217593 


88.3 


24.5 


9.8 


3.3 


228141 


112* 


17.0 


8.5 


6.2 


217678 


43.5 


9.6 


3.8 


4.3 


228230 


24.5 


7.5 


3.8 


3.1 


217689 


119 


27.4 


11.0 


4.1 












217710 


120 


27.7 


11.1 


4.1 












217738 


76.1 


18.0 


7.2 


4.0 












218158 


116* 


18.8 


7.5 


5.8 












218600 


58.8 


11.8 


4.7 


4.7 












HIFI 




492313 


0.88 


0.30 


14.3 


2.8 


1004810 


2.13 


0.50 


4.2 


4.0 


501501 


0.40* 


0.06 


4.0 


6.3 


1007185 


2.50 


0.44 


3.1 


5.3 


504512 


0.14* 


0.05 


3.3 


3.0 


1008602 


2.10 


0.33 


2.5 


6.1 


505339 


0.62 


0.10 


6.3 


6.2 


1008962 


3.16 


0.47 


3.4 


6.3 


507167 


0.72* 


0.07 


4.7 


9.8 


1009245 


1.61 


0.34 


2.4 


4.4 


507673 


0.33 


0.07 


4.3 


4.8 


1009460 


1.85 


0.54 


4.2 


3.2 


515619 


0.25 


0.06 


4.3 


3.6 


1009618 


2.31 


0.62 


5.6 


3.5 


519093 


0.70* 


0.043 


5 


9.1 


1009728 


2.78 


0.51 


4.3 


5.1 


528091 


0.33 


0.06 


4.1 


5.0 


1009804 


2.90 


0.53 


4.4 


5.1 


530769 


0.28* 


0.09 


6.0 


2.8 


1009850 


2.74 


0.58 


4.8 


4.5 


531303 


0.36* 


0.062 


4.1 


5.4 


1009881 


3.20 


0.50 


4.2 


6.0 


531386 


0.65 


0.09 


6.0 


6.8 


1017993 


4.00 


0.71 


7.6 


5.2 


534943 


0.63* 


0.049 


3.3 


12 


1024147 


2.78 


0.60 


6.0 


4.3 


536632 


0.14 


0.06 


3.7 


2.4 


1040932 


1.86 


0.52 


3.7 


3.3 


539134 


0.36 


0.08 


5.3 


4.2 


1045510 


2.30 


0.54 


6.0 


4.0 


556311 


0.44 


0.10 


6.7 


4.1 


1047493 


2.00 


0.47 


3.6 


4.1 


602910 


0.53 


0.10 


3.1 


4.8 


1047982 


3.34 


0.62 


4.1 


5.0 


649022 


0.51 


0.12 


4.0 


4.0 


1052768 


4.50* 


0.75 


3.0 


5.7 


649093 


0.47 


0.12 


4.0 


3.7 


1065139 


3.90 


0.60 


4.6 


6.2 


650170 


0.74 


0.13 


3.8 


5.4 


1066454 


2.00 


0.49 


4.9 


3.8 


692016 


0.39 


0.10 


2.5 


3.7 


1069053 


1.45 


0.30 


3.0 


4.5 


693184 


0.60 


0.13 


3.3 


4.4 


1069207 


1.00 


0.32 


3.2 


2.9 


695857 


0.41 


0.13 


3.3 


2.9 


1069495 


1.35 


0.34 


3.4 


3.7 


697351 


0.95 


0.12 


3.0 


7.7 


1069625 


1.35 


0.32 


3.2 


3.9 


698919 


1.60* 


0.27 


6.8 


5.6 


1069741 


1.48 


0.34 


3.4 


4.1 


706552 


0.55 


0.16 


4.0 


3.2 


1069851 


2.70* 


0.43 


3.9 


5.9 


717678 


0.74 


0.24 


5.3 


2.9 


1076527 


3.80 


0.85 


7.7 


4.2 


741290 


1.24 


0.15 


3.0 


7.9 


1078549 


2.63 


0.51 


4.6 


4.8 


743318 


0.81 


0.15 


3.0 


5.2 


1088059 


2.00 


0.55 


4.6 


3.4 


781479 


1.40 


0.22 


4.4 


6.0 


1088302 


1.65 


0.44 


3.7 


3.4 


785042 


0.85 


0.18 


3.6 


4.4 


1095012 


2.86 


0.69 


6.3 


3.9 


790927 


3.30 


0.52 


5.8 


6.0 


1095149 


1.70 


0.45 


4.1 


3.5 


791560 


1.00 


0.27 


3.4 


3.4 


1098546 


2.15 


0.45 


4.5 


4.4 


795572 


1.60 


0.26 


3.7 


5.7 


1099864 


2.90 


0.62 


6.2 


4.4 


840812 


0.80 


0.21 


4.2 


3.6 


1099903 


1.50 


0.32 


3.2 


4.4 


885659 


1.05 


0.31 


4.4 


3.2 


1108349 


1.96 


0.51 


3.6 


3.6 


972428 


1.70 


0.40 


3.6 


4.0 


1112291 


4.00 


0.48 


3.7 


7.7 


1000929 


3.90 


0.67 


7.4 


5.4 













Notes. Given are the center frequency in the velocity rest frame of SMA2, the line area, the peak intensity, the corresponding signal-to-noise ratio, 
and the line width (FWHM). (6) : blended line. 



are not listed here, but can be checked by comparing the bands 
in H and V polarization or at the edges where the bands overlap. 
In this way, 5 false lines (called ghosts) at 720 GHz are noticed 
which are present in 2a, but not 2b. They arose presumably from 



a feedback in the LO and are no longer present after calibration 
with a newer version of HIPE. Notable are patterns like branches 
where several lines (5-10) appear at once in a narrow frequency 
range (1 GHz), namely at 510, 532, 537, 602, 641, 695, 720, 
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Table 5: Integrated spectral line intensities of all detected chem- 
ical species of the HIFI spectrum for NGC 63341 



Species 


j'Tdv 


Lines 


Species 


j'Tdv 


Lines 




(K-kms- 1 ) 






(Kkms 1 ) 




CH 3 OH 


6016 


2063 


CH 3 CN 


50 


48 


CO 


4715 


31 


HC1 


40 


3 


CH 3 OCH 3 


1486 


637 


N 2 H + 


23 


4 


H 2 


1257 


31 


OCS 


22 


16 


CH 3 OCHO 


582 


558 


H 2 + 


16 


1 


H 2 CO 


428 


105 


NS 


16 


12 


CS 


280 


45 


CN 


15 


6 


C + 


262 


1 


HNC 


15 


4 


C 2 H 5 OH 


227 


134 


CH 2 NH 


14 


17 


NH 3 


226 


13 


CH 


14 


1 


so 2 


207 


109 


OH + 


12 


2 


HCN 


205 


27 


NO 


11 


8 


NH 2 


181 


6 


h 2 ci + 


11 


2 


H 2 S 


168 


32 


C 2 H 


1 1 


6 


HCO + 


152 


10 


HNCO 


10 


25 


HF 


146 


1 


H 2 CCO 


4 


9 


H 2 CS 


141 


78 


HCS + 


4 


4 


SO 


119 


51 


NH 2 CHO 


2 


8 


C 


115 


2 


HC 3 N 


2 


4 


CH + 


108 


3 


SiO 


2 


1 


NH 


103 


3 


SH + 


1 


1 


C 2 H 4 


67 


94 








Sum 


17485 


4217 





Notes. Included are all transitions from isotopologues (except CH 3 OD 
and CH 3 8 OH) and torsionally or vibrationally excited states. 



1009 and 1069 GHz. Other clear patterns like a recurring differ- 
ence between lines could not be discovered. If there are several 
candidates for an unidentified line, the simulation with XCLASS 
helps to constrain them. Any anticoincidence, which means that 
a simulated line in other quantum numbers in other frequency 
bands does not match with the observed spectrum, is an indi- 
cation for a false identification. That is why complete line sur- 
veys are better suited for confident line identifications than small 
windows or single line searches. Another constraint from inter- 
ferometric data is that all lines of a species originate from the 
same source. Since weeds have the highest number of lines, it is 
most probable that the unidentified lines originate from them, or 
one of their isotopologues not yet included in the database. Two 
papers were checked with reported line frequencies of deuter- 
ated meth anol CH 3 OE> (|Walsh et al.|2000| l and the isotopologue 
CH' 8 OH ( Fisher et al. 2007 1. 8 lines could be assigned to deuter- 
ated methanol and 128, especially most of the former reported 
branches, to CH 3 8 OH. Subtracting them, the remaining number 
of unidentified lines is 75 lines (including the 1009 GHz and 
1069 GHz branches) or about 2% of the total of 4428 lines. In the 
SMA bands, 144 transitions are counted of which 27 are uniden- 
tified. That makes 19% and is much more than in HIFI frequency 
range. It is probable that they originate from organic molecules 
in higher excited states. All unidentified lines in the SMA spec- 
tra are detected in both cores, but it is remarkable that, with the 
exception of two, they are all stronger in SMA2 than SMA1. 

5. Discussion 

5. 1 . General remarks 



Schilke et al.] (|2006") analyzed spectra of NGC 63341 with APEX 
at 460 and 810 GHz with a bandwidth of 2.5 GHz. 12 molecules 
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Fig. 4: Derived excitation temperatures in dependence of the 
source size for molecules associated with NGC 63341. 



were identified and roughly fitted, of which 6 show an accept- 
able agreement (same order of magnitude in N) with our results: 



N 2 H + 



and SO.. The 



CH 3 CN, CH 3 OCH 3 , CH3OCHO, HNCO, 
deviation for other molecules can be explained by an over- or un- 
derestimation of the source size and by setting the temperature 
in steps of 50 K. In general, there is a good agreement between 
the HIFI and SMA column densities. The temperatures given in 



Table A. 1 are not well determined, because for many molecules 
only one transition is present. 

Because no error bars are given, it should be pointed out that 
the given parameters for some molecules are less reliable than 
for others. The spectra of strong lines (CO, N 2 H + , HCN, C + ...) 
could be affected by emission in the off-beam position. In some 
cases, the assumption of LTE does not hold, e.g. for HCN, or 
the ambiguity leads to an uncertainty. Reliable values are these 
which are checked consistently with the SMA and have isotopo- 
logues, so that better constraints are available for and N. 

In Fig. |4] the excitation temperature is plotted as a function 
of the source size for all molecules associated with NGC 63341. 
The extended source sizes (9 » beam size) are not included. The 
data points can be roughly fitted by a function fix) — a ■ x b + c 
with an exponent of b = - 1.3 ±0.5 and a constant c — (35±11)K 
for 9 >6" in the outer parts or envelope where the temperature 
gradient becomes flat. This is consistent with models of hot cores 
which predict a steep increase inwards (Rolffs et al.|20lT i. The 
profile could be affected by the fit function, where in the optical 
thick limit, T and 9 are inversely proportional to each other with 
T cc 1/0, but this only concerns few molecules. 

The fractional abundances •X=iVspecies/A / H are distributed 
widely, ranging over eight orders of magnitude from the nitro- 
gen hydrides (X=10~ 12 ) to CO (X=10~ 5 , derived from the iso- 
topologues). The hydrogen c olum n density for NGC 63341 is 



A^h = 3 x 10 24 cm 2 , see Sect. 3.1 The most abundant molecules 



with the highest column densities after CO are CH 3 OH and 
CH3OCH3 (both X=10~ 6 ), HCN and H 2 (both X=10~ 7 ). The 
hydrogen column density can be estimated from the derived 
CO column density of N - 3.7 x 10 19 crrT 2 by using a stan- 
dard fractional abundance of [CO]/[H]=5 x 10 . The result is 
Nu = 7.4 x 10 23 for a source size of 9 = 32" and is lower com- 
pared to the value given above where 9 = 10" is assumed. The 
following differences and similarities are evident when compar- 
ing NGC 63341 with an example of a low-mass protostar IRAS 
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Table 7: Comparison of the fractional abundances in NGC 63341 
with IRAS 16293-2422 and the Hot Core (HQ and Compact 
Ridge (CR) of Orion KL 



Fig. B.l 



Results for most molecules are discussed briefly below, see 
and B.2 for their spectra. They are grouped by their mo- 



Molecule 


NGC 


IRAS 


Orion KL 


Orion KL 




63341 


16293-2422 c 


HC C 


CR e 


CO 


1.2(-5) 


4.0(-5) rf 


9.0(-5) 


1.1 (-4) 


OLOH 


4.7(-6) 


3.0(-7) 


1.4(-7) 


4.0(-7) 


CH 3 OCH, 


1.0(-6) 


2.4(-7) 


8.0(-9) 


1.9(-8) 


HCN 


8.7(-7) 


_ 


3.0(-7) 




H 2 


6.7(-7) 


_ 


5.3(-6) 


1.0(-4) 


OCS 


4.0(-7) 


2.5(-7) 


l.l(-8) 


3.0(-8) 


CH 3 OCHO 


2.4(-7) 


2.0(-7) 


1.4(-8) 


3.0(-8) 


CS 


1.2(-7) 


3.0(-9) rf 


6.0(-9) 


1.0(-8) 


SO 


8.3(-8) 


2.5(-7) 


1.9(-7) 


3.0(-7) 


H 2 S 


8.3(-8) 


9.0(-8) 


2.7C-6/ 




H 2 CO 


5.0(-8) 


6.0(-8) 


7.0(-9) 


4.0(-8) 


so 2 


3.3(-8) 


1.0(-7) 


1.2(-7) 


1.6(-7) 


HCOOH 


1.8(-8) 


6.2(-8) 


8.0(-10) 


1.4(-9) 


C 2 H 5 CN 


1.4(-8) 


1.2(-8) 


3.0(-9) 


5.0(-9) 


H 2 CS 


1.0(-8) 


5.5(-9) 


8.0(-10) 


1.2(-9) 


CH,CN 


9.3(-9) 


1.0(-8) 


4.0(-9) 


5.0(-9) 


HNCO 


5.0(-9) 


9.0(-9) 


5.8(-9) 




C 2 H 


1.4(-9) 


2.2(-10) rf 


5.5(-10) 


5.3(-9) 


HC 3 N 


1.4(-9) 


1.0(-9) 


1.8(-9) 


6.0(-9) 


SiO 


1.0(-9) 


4.5(-9) 


6.0(-9) 


8.0(-9) 


HCO + 


4.3(-10) 


1.4(-9) rf 


1.1 (-9) 


1.0(-9) 


CN 


3.71-10) 


8.0(-ll) rf 


8.0(-10) 


8.0(-10) 



Notes. Notation a(b) indicates a x 10~ b . (c) References: iCazaux et al. 
|2003[ >, ( |Sch6ier et al.|2002]l. (d) referr ing to the cooler, less dense en- 
velope. (eT^eTereTicesT^Sutto^et al. 19951, (Caselli et al. 1993} and 
references therein, (f) from ( |Persson et al.|2007| >. 



16293-2422 and the high-mass star-forming region Orion KL, 
Table The methanol abundance X=5 ■ 1(T 6 in NGC 63341 
is higher than in the Orion HC and IRAS 16293, which have 
the same ratio. As is shown by time-dependent chemical mod- 
els of hot cores ( |Nomur a & Millar 2004), the parent molecule 
methanol (formed on grain surfaces) is stable for timescales 
of at least 10 4 yr. The other organic molecules, dimethyl ether 
and methyl formate, are more abundant in NGC 63341 than in 
Orion, reflecting perhaps that NGC 63341 is more time evolved. 
These second generation molecules start to form in the gas 
phase between evaporated molecules and peak up in abun- 
dance at 10 4 ~ 5 yr. After that, the parent molecules are even- 
tually destroyed. For sulfur-bearing molecules, an anticorrela- 
tion between NGC 63341 and Orion is observable: while H 2 S 



is two orders of magnitude less abundant in NGC 63341, OCS 
is much more abundant than in Orion. The higher depletion of 
H 2 S could be explained by the fact that this molecule is effi- 
ciently destroyed by H 3 + in hot cores and subsequently trans- 
formed into SO and S0 2 fWakelam et~aT][2004) i. Overall, the 
fractional abundance of the sulfur- bearing molecules and their 
sequence in NGC 63341 is more comparable to IRAS 16293 than 
to Orion (except for CS). Some nitrogen-bearing molecules like 
nitriles have about the same fractional abundances in NGC 63341 
as in Orion and IRAS 16293: CH 3 CN with X=l ■ 10" 8 , CN 
with X=4 ■ 10- 10 , HC 3 N with X=2 • 10~ 9 and C 2 H 5 CN with 
X=l • 10~ 8 . When comparing fractional abundances in the lit- 
erature, it should be pointed out that possible uncertainties are 
the adopted H 2 column density and the source size. It needs to 
be further investigated if all these deviations are due to different 
precursor abundances, varying gas temperatures and densities or 
are due to age effects solely. 



ment of inertia, ranging from the lightest diatomic to the heaviest 
complex organic molecules. The electronic state for molecules 
is denoted by 25+1 An, where S is the total electronic spin, A the 
projection of the total electronic orbital angular momenta along 
the internuclear axis and Q = A + 2, where £ is the projection of 
S onto the axis. 



5.2. Atomic and diatomic species 



Ci and Cn, Fig. 
tral carbon fall 



B.la 



into 

3j 



Two fine structure transitions of neu- 
the submillimeter range, 3 / ) i- 3 f > o at 
492 161 MHz and i P 1 -P\ at 809 342 MHz. Up to five veloc- 
ity components can be observed, with the main emission at -7.3 
kms" 1 . The nondetection of I3 C is used to set an upper limit for 
the abundance. The [Cn] 2 / ) 3/2- 2 ^ > i/2 fine-structure transition at 
1.9 THz is one of the most important cooling lines in the ISM 
and one of the most intense single lines (J TAv -262 K-kms -1 ) 
in the HIFI spectrum besides CO. Three components are found, 
corresponding to the ones in Ci. The large source sizes for all 
components indicate a broad distribution. The dust attenuation 
is neglected in the fit of the Cn emission line because it is as- 
sumed that most of the emission originates from the diffuse outer 
layer. Therefore, the derived values for N and T should be taken 
as lower limits, and we notice that some of the spectral features 
could be affected or caused by contamination in the off-beam. 

HF, Fig . B . 1 s The ground transition 7= 1 -0 of hydrogen fluo- 
rideat 1 232476 MHz is seen in absorption (Emprechtinger et al. 
2012 1. Three components at -5, 0.3 and 7 km s" 1 are visible. The 
broad velocity width of the -5 kms 1 component indicates that 
HF is widely distributed and extended over NGC 63341. HF is 
considered to be a good tracer of H 2 ( Neufeld et al.|2005) l, since 
the formation is directly coupled through the exothermic reac- 
tion F + H 2 — > HF + H, and the ratio in the diffuse medium is 
X = 3.6 ■ 10~ 8 . This would correspond to N H = 3 x 10 21 cirT 2 in 
this diffuse layer. 



SO, Fig. B.2d Sulfur monoxide is a radical with two un- 
paired electrons (ground state 3 Z) and a spin-multiplicity of 
3, so that each total angular momentum level J = N + S is 
split into 3 levels according to the projection of S . The nota- 
tion is J,N — > 7 — l,N — 1, where N is the rotational mo- 
menta. Beginning from 7=12, the transitions are strong and can 
be observed up to 7=24. Additionally, the isotopologue 33 SO 
is present for 7=12 and 34 SO for 7=12 up to 7=16. Weaker SO 
transitions are present 43— 12 and 54-23 in Band la and 77-67 up 
to 12]2 - 1 1 12- The emission comes predominantly from SMA2, 
as revealed by the velocity of -7.4 kms -1 and the SMA map. 
The abundance is derived from the isotopologues, but it was nec- 
essary to implement a second component which is optically thin. 

SiO, Fig. |B. lr| Silicon monoxide is viewed as a good tracer 
for outflows and the 7=5-4 transition was observed by Gibb 



et al.| ( |2007[ >. NGC 63341 shows a bipolar outflow with a red 
lobe in the north-east and a blue lobe in the south-west direc- 
tion. With better resolution and additional transitions in the HIFI 
survey the abundance and spatial position can be given more se- 
curely. The 7=12 and 7=13 transitions in the HIFI data are quite 
weak. Therefore the parameters were only manually determined. 
Notable in the SMA data is the narrow - 17 km s _1 component of 
the blue lobe and the wide, more diffuse outflow in the red wing. 
Note that the massive outflows are better traced by CO maps of 
other groups. [Leurini et aL] (p006) analyzed the ratios of several 
CO transitions and derived 1-3 M Q per lobe (dependent on the 
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inclination angle), a dynamical time scale (lobe length divided 
by v max ) of 2 - 4 x 10 3 yr and a total mass entrainment rate of 
M lot = 2 - 15 x l(r 3 M /yr. 



NH, Fig. B.ld The ground state transitions 7=1-0 of the 
radical Imidogen with a spin multiplicity of 3 ( 3 2) and its hyper- 
fine structure are observed. The coupling of the nuclear spin 1=1 
of 14 N with 7 leads to the splitting into F = J + I components. 
To fit it properly, a second broad component had to be included. 
The 7 = I2 - 0] transition is blended by the 7=1 1-10 transition 
ofHCN. 



SH + , Fig. B. lv Sulfoniumylidene has been detected only re- 
cently in Sgr B2"(M) ( |Menten et al.|2011) . The l 2 -0i transition 
and its hyperfine structure are observed tentatively above 3cr in 
absorption and is heavily blended with methanol. The velocity 
matches that of SMA1. A comparative study of SH + absorp- 
tion lines towards other high-mass star-forming regions has been 
made by |Godard etaL] ( |2OT2). 

CH and CH + , Fig. |B.lc| and |B.lb| Th e hyperfine transi- 



tions of methylidine (CH) were analyzed by van der Wiel et al. 
(2010). We included the blended lines from methanol and find 
three components: one broad emission from NGC 63341 and 
two absorption lines from the foreground. The values are in very 
good agreement with their model I. The two lowest transitions 
of methylidynium (CH + ) are observed including the isotope 13 C 
(Lis et al. 2012, in prep.). The three components are approxi- 
mately the same as in CH and show the same abundance. 

NO, Fig. |B.lq| The radical nitrogen monoxide with its dou- 
blet structure ( Z U) is detected in the transitions 2 Ii\ /2, 7 =5.5^1.5 
(e and f) up to 2 Yl\/2, 7 =8.5-7.5 (e and f). The hyperfine struc- 
ture is not resolved. 



NS, Fig. B.2b The electronic ground state of nitrogen sul- 
fide is 2 n, which splits the rotation levels due to the unpaired 
electron. All transitions from from 7=11 to 19 are detected. 
Noticeable is that the fit for the H1/2 states is good, but the TI3/2 
states are underexcited. The optical depth is small, so the source 
size was fixed at 9 = 3". 

CN, Fig. |B.ll| Due to the weak emission and small number 
of transitions, there is a large uncertainty of the abundance de- 
pending on the velocity width. A width of 4 km s _I would result 
in an abundance higher by a factor of 2. Here the best fit is taken. 
The 7 = 203 - I0 2 transition at 226.8 GHz in the SMA data is 
seen in absorption and resembles well the dust emission. The as- 
signment to SMA1 due to visr = -6.4 kms _I is confirmed by 
the map. 

CO and CS, Fig. B.lm and B.2a Carbon monoxide is the 
most abundant molecule after H 2 in the ISM. In the HIFI spec- 
trum, CO has the most intense lines, which are optically thick. 
The line profile shows strong outflows with widths of over 
80 kms 1 and self-absorption. The same is true for carbon 
monosulfide, so no direct abundance determination is possible. 
Therefore only the four detected isotopologues of CS and CO are 
considered, with a focus on C ls O and C n O by Ceccarelli et al. 
(2012), in prep. For CS, a multi-component approach was tried 
with two components from the hot cores and one extended emis- 
sion. But because the first and third component have the same 
offset velocity, there might be a redundancy. 



and a narrow line width had to be taken. Russ eil et al.| ([2010) 
find two components also by observing the hyperfine structure 
of the loi - O12 line, at -7.3 and -10 kms -1 , but with large dif- 
ferences between them concerning T (127 K to 5 K) and Ay (4 
to 19 km s _1 ). Another possibility could be self-absorption. 



HNC and HCN, Fig. B.lo and B.ln Hydrogen isocyanide is 



an isomer of hydrogen cyanide (HCN). Transitions from 7=6 to 
10 are present for HNC and 7=6 to 12 for HCN. As for N 2 H\ 
a two component model was taken. The line profile is chang- 
ing and affected by self absorption. HCN has a high critical 
density (n cr > 10 9 for 7=9-8 at 797 GHz) and shows an out- 
flow in the red wing of the asymmetric line profile with self- 
absorption, which indicates an infall. Also vibrationally excited 
lines (v 2 = 1,2) of HCN are observed. Both aspects, infall and 
hot HCN, are characteristic for evolved cluster-forming regions 
(Rolffs et al. 201 1 1. A uniform fit with the isotopologues of HCN 
is included for the completeness of the survey. For HNC, no 
other isotopologues are found apart from the HN 13 C 7=6-5 line. 



HCO + , Fig. B.2p The rotational transitions of the formyl ion 
are relatively intense and present up to band 5a (Ceccarelli et al. 
2012, in prep.). At first, one component including H 13 CO + was 
tried, but the fitting analysis revealed that either the lower-7 or 
the higher-7 transitions are consistent but not all together. This 
is an indication that either LTE is not satisfied, or several com- 
ponents with different excitation conditions overlap. Assuming 
two components with different optical depths, the fit could be 
improved. At higher transitions the component with the higher 
optical depth dominates. 



OCS, Fig. B.2c Carbonyl sulfide is detected in the rotational 
transitions from 7=40 to 53. Because at first the source size was 
unknown and no isotopologues are detected in the HIFI spec- 
trum, the abundance was underestimated due to a low opacity 
(t < 0.3). With the detection of three isotopologues in the SMA 
bands the column density was corrected upwards. The emission 
of 13 CS comes predominantly from SMA1. 



HC 3 N, Fig. B.2t Cyanoacetylene lines are weak in the HIFI 



data, but more intense in the SMA data. The fit for SMA1 gives 
a large line width of 9 km s _1 , which is taken for the HIFI model. 

5.4. Symmetric top molecules 



CH 3 CN, Fig. B.2m| Methyl cyanide is considered to be a good 
tracer of hot cores ( |Olmi et al.|1993[ ). Rotational transitions from 
7=27 to 39 with /T-components K=0-6 are detected in the HIFI 
spectrum. |Thorwirth et al.| ([2003) fitted the 7=12-11 transition 
with an unknown source size resulting in 7=145 K and N=5.2 x 
10 15 cm -2 . The temperature is in good agreement with our value 
of 154 K. Additionally, weak (< 3cr) bending modes vs=l are 
observed in band la. The small source size and high temperature 
confirms that methyl cyanide probes the central part of hot cores. 



NH 3 , Fig. B.lf The ground transitions of ammonia are de- 



tected mostly in absorption. The overlapping emission and ab- 
sorption components could not be fitted consistently, and the 
column density for the emission component associated with 
NGC 63341 is therefore probably underestimated. 



5.3. Polyatomic linear molecules 



N 2 H + , Fig. B.lg Four emission lines from 7=6 to 10 of di- 
azenylium are detected which exhibit an asymmetric line pro- 
file (Ceccarelli et al. 2012, in prep.). The intensities could not 
be fitted properly with a single component. The two component 
model shows a much better fit, but the values are very uncertain 



5.5. Inorganic asymmetric molecules 



NH 2 , Fig. B.le Azanyl is an asymmetric top with a hyper- 



fine structure. Four transitions are detected including the ground 
state transition. The emission component can be associated with 
NGC 63341. As in the case of NH, a broad component (Av=30 
kms 1 ) is added to fit the absorption. 
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H 2 C1 + and HC1, 



Fig. 



B.lu and B.lt 



chloronium were analyzed by Lis et a 



Hydrogen chloride and 
" ((20T0b. The 7=2-1 at 



CH3OH, Fig. B.2g and B.2h Methanol is the strongest weed 



1251 GHz in band 5b of HC1 is additionally detected, showing 
a second component in absorption from the colder and more ex- 
tended envelope at -10 kms~' with a broad width comparable 
to H 2 C1 + . The emission component can be assigned to the cen- 
ter of NGC 63341. The abundance for this component differs by 
50% compared to ( Lis et al.|2010| l because the source size is not 
known and the lines are optically thick and affected by the ab- 
sorption feature. For H 2 C1 + , the blended lines of dimethylether 
are included in the fit (blue line). Its velocity and line width leads 
to the conclusion that it originates from a diffuse foreground 
molecular cloud. 



H 2 + 



Fig. B.lj 



Oxidaniumyl was analyzed by Ossenkopf 
et al. (2010 1. Additionally, another weak fine structure transition 
was observed at 1140 GHz from the 7=0.5 level. The derived 
abundance is above the lower limit given by |Ossenkopf et al. 



( 2010| l. Its radial velocity and line width leads to the conclusion 
that it is not associated with NGC 63341. 



FLO, Fig. B.li Water was analyzed by Emprechtinger et al. 



(2010 2012) . It has complex line shapes with many components 



like optically thick emission and broad outflows, therefore no fit 
was attempted. Instead, the optically thin isotopes were included 
to get an estimation of the abundance. Most fits including deuter- 
ated water can not reproduce the line strength, so that a non-LTE 
calculation method is necessary to model them properly. 

Hydrogen sulfide shows 3 components. The 
-6.1 kms -1 is affected by two absorption lines, 



H 2 S, Fig. 



B.lw 



emission one at 

of which one is from the outflow. While the emission peak could 
be constrained through the isotopologues, the absorptions fea- 
tures could not be reproduced well. 



5.6. Organic asymmetric molecules 



H 2 CO. Another less efficient 



H 2 CO + H. 



Besides 



H 2 CO, Fig. B.21 Formaldehyde is a nearly prolate top and the 
asymmetry splits every K + level into doublets. For astro- 
physical applications, due to its structure this molecule is a good 
probe of the kinetic temperature and volume density and traces 
the inner and outer part of hot cores ( Mangum & Woot ten|1993| l. 
Interstellar formaldehyde is formed by hydrogenation of solid 
CO: H + CO — > HCO + H 
pathway in the gas phase is CH 3 + O 
of the weeds, this organic molecule has the most intense lines 
throughout the HIFI spectrum from band la to 4b. The rotational 
transitions range from 7=7 to 15 with K-Q to K-6, and the 
isotopologues H 2 3 CO and H2C ls O are observed. The fits with 
XCLASS show two problems: one component is optically thick, 
and all the K components could not be described by a single ro- 
tation temperature. The K-Q, 1 do not match with K-2-6 com- 
ponents consistently. From the SMA data the 3o - 2o transition 
shows that the K-Q value is lower in intensity than the K-2 com- 
ponent. This hints to self absorption in SMA1 or optical thick- 
ness. 



H 2 CCO and H 2 CS, Fig. [B~2"o] and [R2f| Ethenone and thio- 
formaldehyde are nearly prolate tops with a spectrum similar to 
H 2 CO. For H 2 CCO, only the K=\ components from 7=24 to 27 
are detected, the other K-Q and K-3 components are included to 
set an upper limit. Since the transitions are weak and the optical 
depth is small, and Av needed to be fixed. For H 2 CS rang- 
ing from 7=14 to 23 with K=Q to K=5 transitions, a good fit is 
obtained where all K components can be described consistently 
with one excitation temperature. 



in the HIFI spectrum and accounts for about half of all detected 
lines. One component is optically thick, which can be derived 
from the fit of the 13 CH 3 OH isotopologue, and the other compo- 
nent accounts for the broad emission. The model fits the observa- 
tion quite well, except of the most intense lines. As can be seen 
in the SMA spectra, the XCLASS LTE model cannot reproduce 
the methanol lines due to a high optical depth. Furthermore, high 
energy lines > 500 K are underestimated, so that a 3rd compo- 
nent might be necessary. 



NH 2 CHO, Fig. B.2s Formamide is only detected weakly in 



band la. The optical depth is small, and the source size is deter- 
mined from the SMA map which shows that the emission comes 
mainly from SMA2. 



C 2 H 4 0, Fig. B.2n Ethylene oxide is surprisingly only de- 
tected in SMA2. As for formamide, many weak lines are ob- 
served. The source velocity v L s R in the HIFI data is consistent 
with the general source velocity of SMA2. 



HNCO, Fig. B.2r Isocyanic acid is almost linear, but asym- 
metric and shows relatively weak emission. The very high tem- 
perature and small source size may trace the hot core SMA2. 

CH3OCH3, Fig. |B.2i| Dimethyl ether has two methyl 
groups, and similarly to methanol, a splitting into 4 sublevels 
(AA,EE,AE,EA). It is, after methanol, the organic molecule with 
the most dense and intense lines in the entire HIFI spectrum. 
The laboratory measured spectrum from the 13 CH 3 OCH, iso- 



topologue are taken from (Koerber et al. 2010), but are not yet 
publicly available. The SMA spectra show that the dominant 
emission arises from SMA1 and that they are weak associations 
with the outflows. 



CH3OCHO, Fig. |B.2j Methyl formate has the densest spec- 



(Brouillet et al. 2012, in prep.) and the highest number of lines 
in both SMA bands. The SMA maps show an association with 
the outflows. 



6. Conclusions 

In this article, observational data from the HIFI instru- 
ment aboard the Herschel Space Observatory and from the 
Submillimeter Array are analyzed for the high-mass star- 
forming region NGC 63341. We employ a LTE radiative transfer 
model in order to get insights into the physical structure of the 
source. The molecular spectral line survey gives an overview of 
the chemical inventory of the molecular species and their excita- 
tion, whereas the high-resolution spectrometer allows tracing the 
kinematics of the different velocity components. The interferom- 
eter maps reveal the source morphology and help to distinguish 
the embedded cores and resolve their sizes. The results fit well 
with the overall picture of a high-mass young stellar object and 
hot molecular cores described by |van Dishoeck & Blake ( 1998 1. 

Two main software tools, XCLASS and MAGIX, are used 
to demonstrate that modeling and handling of the vast HIFI 
data sets is possible. These tools allow to decompose the spec- 
trum into contributions of individual molecules and their iso- 
topologues. By assuming LTE condition and employing the LM 
fitting algorithm, a fast extraction of information is possible. 
In order to take line blends into account, a model for every 
molecule is produced, which is consistent over multiple HIFI 
bands. It should be clear that for other excitation conditions than 
LTE and for optically thick lines like CO this method reaches 
its limitations. Sophisticated 3-d radiative transfer codes such as 
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RADMC-3E0and LIME <|Brinch & Hogerheijde|2010) with the 
possibility of non-LTE calculations will be used in the future in 
order to model NGC 63341 into more detail, taking into account 
the geometry and non-homogeneous distribution of tempera- 
ture and abundances. With the improved capabilities of ALMA 
(Atacama Large Millimeter Array), it will be possible to study 
the molecular chemistry and kinematics on sub-arcsecond scales 
(< 1000 AU) and to follow the fragmentation of protostellar 
cores with much improved mass sensitivity. 

The molecular survey revealed a variety of chemical species, 
including saturated organic molecules. The astrochemical ques- 
tion is: how can these abundances be explained and how are 
they related to each other? The next step is therefore a careful 
comparison of the column densities presented here with predic- 
tions of astrochemical models and detailed comparisons with the 
molecular abundances and abundance ratios in other CHESS tar- 
gets. 
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Appendix A: SMA spectra and tables 



Table A.l: XCLASS fit results from SMA for NGC 63341 



Source 


Species 





T„ 




N 


Av 


Vlsr 


Notes 






(") 


(K) 


(cm- 2 ) 


(km/s) 


(km/s) 




SMA1 


CN 


2 


3 


0.2 x 


10 ls 


2.2 


1.3 




SMA1 


CN 


2 


3 


0.8 x 


10 15 


6.2 


-4.6 




SMA2 


CN 


2 


3 


0.7 x 


10 15 


6.9 


-0.6 




SMA2 


CN 


2 


3 


1.5 x 


10' 5 


3.2 


-4.5 




SMA1 


DCN 


3.5 


359 


1.6 x 


10 15 


5.8 


-4.9 


2 


SMA2 


DCN 


5.0 


146 


0.6 x 


10 15 


6 


-7.8 


2 


SMA1 


SiO 


1.7 


40 


l.Ox 


10 15 


8 


-3.2 


2 


SMA2 


SiO 


2.0 


40 


1.0 x 


10' 5 


6 


-5.2 


2 


SMA1 


H 2 S 


3 


62 


1.4 x 


10 17 


5 


-4 


2,3 


SMA2 


H 2 S 


3 


147 


0.4 x 


10 17 


6.5 


-6 


2,3 


SMA1 


CS 


2.4 


40 


1.6 x 


10 18 


5.1 


-3.2 


2,3,5 


SMA2 


CS 


2.2 


32 


6.6 x 


10 17 


2.7 


-6.4 


2,3,5 


SMA1 


OCS 


2.5 


137 


2.3 x 


10 18 


4 


-5.6 


1 


SMA2 


OCS 


2.5 


126 


1.1 X 


10 18 


4 


-7 


1 


SMA1 


CH 3 OH 


4 


78 


2.4 x 


10 19 


3.6 


-5.2 




SMA2 


CH3OH 


4 


104 


2.0 x 


10" 


3.0 


-7.5 




SMA1 


CH^OCHj 


3.0 


78 


2.3 x 


10 18 


4.9 


-5.6 




SMA2 


CH~OCH 3 


3.0 


98 


0.7 x 


10 18 


5.4 


-7.7 




SMA1 


CH^OCHO 


2.3 


104 


7.7 x 


10 17 


4.6 


-6.1 




SMA2 


CH^OCHO 


3 


115 


9.6 x 


10" 


3.8 


-7.6 




SMA1 


H 2 CO 


3.0 


104 


1.6 x 


10 17 


5.7 


-5.5 




SMA2 


H 2 CO 


3.0 


138 


1.3 x 


10 17 


4.6 


-7.5 




SMA1 


NH 9 CHO 


2 


25 


4.7 x 


10 15 


5.2 


-4.7 


2,3 


SMA2 


NH 2 CHO 


2 


45 


5.0 x 


10 15 


4.0 


-6.9 


2,3 


SMA1 


HC 3 N 


2.8 


208 


4.2 x 


10 15 


9.4 


-5.7 




SMA2 


HC^N 


2.4 


223 


1.4 x 


10 15 


5.8 


-6.6 




SMA1 


C 2 H 3 CN 


2.0 


147 


5.2 x 


10 15 


6.0 


-6.3 




SMA1 


C 2 H 5 CN 


2.2 


280 


4.2 x 


10 16 


6.6 


-5.5 




SMA2 


C 2 H 5 CN 


1.0 


123 


1.6 x 


10 16 


4.1 


-6.7 




SMA1 


CH^CHjCO 


2.6 


73 


2.1 x 


10 16 


5.6 


-5.4 




SMA2 


CH3CH3CO 


2.6 


159 


2.9 x 


10 16 


4.8 


-7.4 




SMA2 


HCOOH 


2.0 


258 


5.5 x 


10 16 


4.3 


-9 


2 


SMA2 


(CH 2 OH) 2 


2.2 


217 


2.6 x 


10 16 


3.4 


-8.2 




SMA2 


C 2 H 5 OH 


2.5 


77 


1.4 x 


10 17 


4.9 


-7.9 




SMA2 


C 2 H 4 


2.1 


335 


6.7 x 


10 16 


5 


-6.8 




SMA2 


SO 


2.5 


70 


1.7 x 


10 17 


4.2 


-8.7 


2,6 



Notes. (1) including only isotopologues 18 OCS, 13 CS, OC 33 S. (2) 1 or few lines. (3) blended. (4) including 13 CH 3 OH. (5) including only 13 C 34 S. 
(6) including only 33 SO. 
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Table A.2: List of identified lines in the SMA bands of NGC 63341 in order of increasing rest frequency 



Molecule 



Transition 



LSB 

Frequency (MHz) 



Molecule 



Transition 



USB 

Frequency (MHz) 



H 2 S 

C 2 H 4 

C 2 H 4 

18 OCS 

CH3OCHO 

CH3OCHO 

CH3OCHO 

CH3OCHO 

CH 3 CH3CO,v=l 

CH 3 CH 3 CO,v=l 

CH 3 CH3CO,v=l 

CH 3 CH3CO,v=l 

CH3OH 

CH3OCHO 

CH3OCHO 

CH 3 CH3CO,v=0 

CH 3 CH3CO,v=0 

CH 3 CH 3 CO,v=0 

CH 3 CH3CO,v=0 

CH 3 CH 3 CO,v=0 

CH 3 CH3CO,v=0 

13 CH 3 OH 

CH 3 CH 3 CO,v=0 

CH 3 CH3CO,v=0 

CH3OCHO 

SiO 

CH3OCHO 
CH3OCH3 
CH3OCHO 
DCN 

CH3OCHO 

CH3OH 

CH3OCHO 

13 CH 3 OH 

(CH 2 OH) 2 

(CH 2 OH) 2 

CH3OH 

C 2 H 5 OH 

C 2 H 5 OH 

(CH 2 OH) 2 

C 2 H 4 

C 2 H 4 

CH3OH 

CH 3 CH 3 CO,v=l 

CH 3 CH 3 CO,v=l 

C 2 H 5 OH 

33 SO 

CH3OH 

CH 3 CH3CO,v=l 

CH 3 CH3CO,v=l 

CH 3 CH 3 CO,v=0 

CH 3 CH 3 CO,v=0 

CH3OCHO 

CH 3 CH 3 CO,v=0 

CH 3 CH 3 CO,v=0 

CH 3 CH 3 CO,v=0 

CH 3 CH3CO,v=0 

13 CS 

HjCO 

CH3OCHO 

CH3OCHO 

CH3OCHO 

HC 3 N,v=0 

C 2 H 4 

C 2 H 4 

(CH 2 OH) 2 



C 2 H 3 CN 

C 2 H 3 CN 

C 2 H 3 CN 

C 2 H 3 CN 

CH 3 CH 3 CO,v=l 

CH 3 CH3CO,v=l 

CH3OH 

C 2 H 3 CN 



NH 2 CHO 

C 2 H 5 OH 

H 2 CO 

CH3OCH3 

CH3OCH3 

C 2 H 3 CN 



C 2 H 3 CN 
C,H,CN 



220-2 11 

165 11 - 164 12 

166 11 - 165 12 
19- 18 

32 6 27 1 - 32 5 28 1 
32 6 27 - 32 5 28 
18 2 16 2 - 17 2 15 2 

18 2 160- 17 2 15 
194 160- 184 15 1 

19 3 160- 184 15 1 
194 160- 183 15 1 
19 3 160- 18 3 15 1 
5 100-4200 

17 3 14 5 - 16 3 13 5 
200202- 19 1 19 1 
19 4 16 1 - 18 3 15 1 
19 3 16 1 - 18 4 15 1 
194 16 1 - 183 152 
19 3 16 1 - 18 4 15 2 
194 160- 18 3 15 1 
19 3 160- 184 15 1 

14 1 13 0- 13 2 12 
194 160- 183 150 

19 3 160- 184 15 
30 4 26 - 30 4 27 
5-4 

35 9 27 1 - 35 8 28 1 
224 19 3 - 223 203 
32 9 24 - 32 8 25 
32-22 

30 4 26 - 30 3 27 
6 1-1 1-72-11 
174 133 - 164 123 
10280-9370 

24 1 24 - 23 1 23 1 

24 24 - 23 23 1 

16 1 00 - 15 3 00 
5 1 4 1 -4040 

25 3 22 - 24 4 20 1 
21 2 19 1 - 20 2 18 
154 11 - 153 12 

15 5 11 - 154 12 
15 6 1 1 - 165 1 1 
202 180- 192 170 

20 3 18 0- 19 3 17 
5332-4222 
568-457 

20 1 - 20 
20 2 18 0- 19 3 17 1 
20 3 18 0- 19 2 17 1 
20 3 18 1 - 19 2 17 2 
20 2 18 1 - 19 3 17 2 
174 135 - 164 125 
20 2 18 0- 19 3 17 1 
20 3 18 0- 19 2 17 1 

20 3 18 0- 19 2 17 
202 180- 193 170 
18- 17 

3 3-202 

31 9 23 - 31 8 24 

17 3 14 2 - 16 3 13 2 
17 3 140- 163 13 
24-23 

143 11 - 142 12 

144 11 - 143 12 
224 190-21 4 18 1 
243 21 - 23 3 20 

23 7 17 - 22 7 16 
23 7 17 - 22 7 16 
23 7 16 - 22 7 15 
23 6 18 - 22 6 17 
23 6 17 - 22 6 16 
11 11 1 0- 10 1000 
11 11 00- 10 10 1 

4 2 0-3100 
23 5 19 - 22 5 18 
23 5 18 - 22 5 17 
10 1 9-9 1 8 
5322-4232 
322-22 1 

23 3 21 1 - 23 222 1 
23321 0-232220 
23 10 13 -22 10 12 
23 10 14-22 10 13 

21 5 16 2-21 4 17 2 
23 4 20 - 22 4 19 

23 3 21 - 223 20 
27721 0-275220 



216710.437 

216741.468 

216741.712 

216753.488 

216758.969 

216809.014 

216830.197 

216838.889 

216881.337 

216881.337 

216881.337 

216881.337 

216945.521 

216958.834 

216962.989 

216974.426 

216974.426 

216974.464 

216974.464 

217022.509 

217022.509 

217044.616 

217070.504 

217070.504 

217077.079 

217104.98 

217165.366 

217189.668 

217215.847 

217238.631 

217266.549 

217299.205 

217312.626 

217399.55 

217449.995 

217450.27 

217525.002 

217548.152 

217549.328 

217587.548 

217615.673 

217615.72 

217642.678 

217744.747 

217744.747 

217803.689 

217832.642 

217886.504 

218023.182 

218023.182 

218091.448 

218091.448 

218108.438 

218127.207 

218127.207 

218162.929 

218162.929 

218198.998 

218222.192 

218259.594 

218280.9 

218297.89 

218324.723 

218335.626 

218335.634 

218371.495 

218389.97 

218398.555 

218398.555 

218398.555 

218402.435 

218402.451 

218402.853 

218404.816 

218440.063 

218451.297 

218452.357 

218459.653 

218461.226 

218475.632 

218491.914 

218494.39 

218519.997 

218519.997 

218554.499 

218573.646 

218585.072 

218593.341 



C,H 5 OH 
CN 

CH 3 CH 3 CO,v=l 

CH 3 CH 3 CO,v=l 

CH3OCHO 

CH3OCHO 

13 CH 3 OH 

CH 3 CH3CO,v=0 

CH 3 CH 3 CO,v=0 

CH 3 CH3CO,v=0 

CH 3 CH 3 CO,v=0 

CH3OCHO 

CH3OCHO 

CH 3 CH 3 CO,v=0 

CH 3 CH 3 CO,v=0 

CH3OCHO 

CH3OCHO 

CN 

CN 

CH 3 CH 3 CO,v=0 

CH 3 CH 3 CO,v=0 

CN 

CN 

CH3OCHO 

CH3OH 

CH3OCHO 

CH3OCHO 

CH3OCHO 

CH 3 CH3CO,v=l 

CH 3 CH3CO,v=l 

CH 3 CH 3 CO,v=l 

CH 3 CH 3 CO,v=l 

CH3OH 

CH3OCHO 

CH3OH 

C,H,OH 

13 C 3 *S 

HC 3 N,v=0 

CH 3 CH3CO,v=l 

CH 3 CH 3 CO,v=l 

CH3OCHO 

(CH 2 OH) 2 

CH3OCHO 

NH 2 CHO 

C 2 H 5 OH 

C 2 H 5 CN,v=l 

CH 3 OCH 3 

CH 3 OCH 3 

C,H 5 CN 

CH3OH 

CH 3 CH3CO,v=l 

CH 3 CH3CO,v=l 

C 2 H 3 CN 

C 2 H 3 CN 

CH 3 CH3CO,v=0 

CH 3 CH3CO,v=0 

CH 3 CH 3 CO,v=0 

CH3CH 3 CO,v=0 

C,H 3 CN 

C 2 H 3 CN 

CH 3 CH3CO,v=0 

CH 3 CH 3 CO,v=0 

CH 3 CH 3 CO,v=0 

CH 3 CH3CO,v=0 

CH3OCHO 

C,H 3 CN 

C 2 H 3 CN 

CH 3 CH 3 CO,v=0 

CH 3 CH 3 CO,v=0 

HC 3 N,v 7 = l 

CH 3 OCHO 

CH 3 OCHO 

C 2 H 5 OH 

C,H 5 OH 

CH 3 OCHO 

C,H 3 CN 

C 2 H 3 CN 

C 2 H 3 CN 

OC 33 S 

C 2 H 3 CN 

CH 3 OCHO 

CH 3 OCHO 

CH 3 OCHO 

HC 3 N,v 7 =l 

CH 3 OCH 3 

CH 3 OCH 3 

CH 3 OH 

C 2 H 5 CN 

HCOOH 



10292-9 1 82 
2 2 2-1011 
20 3 17 0- 19 4 16 1 
20 4 17 0- 19 3 16 1 
20 2 19 1 - 19 2 18 1 
20 2 19 0- 19 2 18 
125 8 0- 13 490 
20 4 17 1 - 19 3 16 1 
20 3 17 1 - 19 4 16 1 
20 4 17 1 - 19 3 16 2 
20 3 17 1 - 19 4 16 2 
20 1 19 2 - 19 1 18 2 
20 1 19 0- 19 1 18 
20 4 17 0- 19 3 16 1 
20 3 17 0- 19 4 16 1 
20 2 19 1 - 19 1 18 2 
20 2 19 0- 19 1 18 
2 3 3-1022 
2034-1023 
20 4 17 0- 19 3 16 
20 3 17 0- 19 4 16 
2 3 2-1022 
2 3 3-1023 
66 14-55 14 
14-200- 13 30 
25 9 17 1 - 25 8 18 1 
19 2 17 2 - 18 2 162 

19 2 170- 18 2 160 

20 3 17 1 - 19 4 16 1 
20 4 17 1 - 19 3 16 1 
20 4 17 1 - 19 3 16 2 

20 3 17 1 - 19 4 16 2 

21 1 00-21 00 

25 9 16 2 - 25 8 18 1 
12-1 00- 11 200 
13 3 10 1 - 12 3 9 1 
5-4 

25-24 

21 3 19 0- 203 18 
21 2 19 0- 202 18 
21 21 2-200202 
224 180-21 4 17 1 
18 3 15 3 - 17 3 14 3 
110 11-10 10 
18 5 13 2 - 18 4 14 2 

26 1 26 1 - 25 1 25 1 
47 11 36 1 -46 12 35 1 
47 11 37 3 -46 12 34 3 
25 3 22- 243 21 

16 1 1 - 152 1 
21 3 19 0- 203 18 1 
21 2 19 0-202 18 1 
24 7 18 - 23 7 17 
24 7 17 - 23 7 16 
21 3 19 1 - 20 3 18 1 
21 2 19 1 -20 2 18 1 
21 3 19 1 - 20 3 18 2 
21 2 19 1 - 20 2 18 2 
24 6 19-23 6 18 
24 6 18 - 23 6 17 
21 2 19 0-202 18 1 
21 2 19 0- 203 18 1 
21 3 19 0- 202 18 1 
21 3 19 0- 203 18 1 
24 9 15 - 24 8 16 
24 5 20 - 23 5 19 

24 5 19 - 23 5 18 

21 3 19 0- 203 18 
21 2 19 0- 202 18 

25 1 -24-1 

31 4 27 2 - 31 3 28 2 
24 9 15 2 - 24 8 16 2 
13 1 12 1 - 12 1 11 1 
13 3 100- 123 90 
31 4 270-31 3280 
24 11 13 -23 11 12 
24 11 14-23 11 13 
24421 -23420 
19-18 

24 4 20 - 23 4 19 
24 9 16 1 - 24 8 17 1 
18 3 15 5 - 17 3 14 5 

24 9 16 - 24 8 17 

25 -1 - 24 1 
263240-25421 

26 3 24 1 -25 4 21 1 
21 000-20300 
25 2 23 - 24 2 22 

10 2 8 - 9 2 7 



226661.701 

226679.311 

226692.48 

226692.48 

226713.06 

226718.688 

226716.392 

226784.639 

226784.639 

226784.67 

226784.67 

226773.13 

226778.786 

226832.055 

226832.055 

226856.872 

226862.257 

226874.191 

226874.781 

226879.393 

226879.393 

226887.42 

226892.128 

226913.736 

226939.41 

226974.039 

227019.55 

227028.121 

227032.073 

227032.073 

227043.496 

227043.496 

227094.747 

227106.67 

227229.511 

227294.752 

227300.506 

227418.905 

227560.862 

227560.862 

227561.741 

227587.074 

227599.261 

227606.176 

227606.079 

227624.735 

227668.331 

227675.733 

227780.972 

227814.528 

227836.273 

227836.273 

227897.607 

227897.608 

227903.922 

227903.922 

227903.953 

227903.953 

227906.683 

227906.709 

227939.374 

227939.374 

227939.374 

227939.374 

227954.452 

227966.032 

227967.589 

227974.758 

227974.758 

227977.277 

227986.255 

227994.545 

227891.911 

228029.05 

228057.899 

228087.245 

228087.245 

228104.614 

228151.943 

228160.305 

228205.832 

228211.291 

228270.535 

228303.174 

228421.737 

228422.878 

228467.902 

228483.136 

228544.168 
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Appendix B: HIFI fitted curves of all molecules 

Examples of selected, fitted spectra for each molecule. In red 
is the model, in blue are mostly annotations or blended lines. If 
only few transitions of a molecule are present, they are all shown 
in the plot. For weeds or molecules with many lines, no quantum 
numbers are given and only selected frequency ranges over dif- 
ferent bands are presented. Heavy blended lines are excluded for 
simplicity. All intensities are given in units of Kelvin. The quan- 
tum numbers refer to upper state - lower state. For absorption 
lines, this sequence should be read reversed. 

The entire observed and simulated spectrum for each 
HIFI band will be made available online in the CLASS 
format on the CHESS website (http://www-laog.obs. 
ujf-grenoble.fr/heberges/hs3f/), as soon as the HIPE 8 
reprocessed data are available. Additionally, the modeled spec- 
trum for each molecule (including isotopologues and vibra- 
tionally excited states) can be loaded. It should be noticed that 
for the simulated spectrum of only one molecule, the blending 
by other molecules is not taken into account. This results in a 
deviation for the fits of absorption lines when a nearby emis- 
sion line contributes to the continuum. In XCLASS, at first all 
the emission components are calculated and then the absorption 
components. 
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